Superoperator Analysis of Entanglement in a Four-Qubit Cluster State 
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In this paper we utilize superoperator formalism to explore the entanglement evolution of four- 
qubit cluster states in a number of decohering environments. A four-qubit cluster state is a resource 
for the performance of an arbitrary single logical qubit rotation via measurement based cluster state 
quantum computation. We are specifically interested in the relationship between entanglement 
evolution and the fidelity with which the arbitrary single logical qubit rotation can be implemented 
in the presence of decoherence as this will have important experimental ramifications. We also 
note the exhibition of entanglement sudden death (ESD) and ask how severely its onset affects the 
utilization of the cluster state as a means of implementing an arbitrary single logical qubit rotation. 

PACS numbers: 03.67.Mn, 03.67.Bg, 03.67.Pp 



I. INTRODUCTION 

Entanglement is a uniquely quantum mechanical phe- 
nomenon in which quantum systems exhibit correlations 
above and beyond what is classically possible. Entangled 
systems are thus an important resource for many quan- 
tum information processing protocols including quantum 
computation, quantum metrology, and quantum commu- 
nication [l[ . In the area of quantum computation, certain 
entangled states play a unique role as the basic resource 
for measurement-based quantum computation. The clus- 
ter state in particular allows for quantum computation 
to proceed via single qubit measurements after creation 
of the cluster state @] . 

An important area of research is to understand the 
possible degradation of entanglement due to decoherence. 
Decoherence, stemming from unwanted interactions be- 
tween the system and environment, is a major challenge 
confronting experimental implementations of quantum 
computation, metrology, and communication [3j. Deco- 
herence may be especially detrimental to highly entan- 
gled states jj] and, indeed, much work has been done on 
studying the effects of decoherence on cluster states Q. 

An extreme manifestation of the detrimental effects of 
decoherence on entangled states is "entanglement sudden 
death" (ESD), in which entanglement within a system is 
completely lost in a finite time [f| 0] despite the fact 
that the loss of system coherence is asymptotic. This 
aspect of entanglement has been well explored in the case 
of bi-partite systems and there are a number of studies 
looking at ESD in multi-partite systems [l|-[l3j including 
the four qubit cluster state flij . In addition, there have 
been several initial experimental ESD studies flo| . 

In this paper we study the entanglement evolution of a 
four-qubit cluster state which can be used as the basic re- 
source to perform an arbitrary single (logical) qubit rota- 
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tion via cluster state quantum computation. We analyze 
the effects of various decoherence models on the entan- 
glement of the pre-measurement state and compare the 
entanglement behavior to the accuracy with which the 
decohered state can be used to implement the desired ar- 
bitrary single qubit rotation. To completely characterize 
the effects of decoherence we make use of superoperator 
representations and aspects of quantum process tomog- 
raphy. Quantum process tomography is an experimental 
protocol which is used to completely determine (open) 
system dynamics. The information gleaned from quan- 
tum process tomography can, in turn, be used to deter- 
mine a wealth of accuracy measures. One would expect 
that the proper working of cluster state based quantum 
computation would be strongly dependent on the amount 
of entanglement present in the pre-measurement cluster 
state. Thus, an explicit analysis of the strength of this 
dependence, especially when attempting to perform basic 
computational gates, is essential for experimental imple- 
mentations of cluster state quantum protocols. 



A secondary aim of this paper is to analyze the ef- 
fect of ESD on the implementation of the single logical 
qubit rotation. The ESD phenomenon is of interest on 
a fundamental level and important for the general study 
of entanglement. However, it is not yet clear what the 
effect of ESD is on quantum information protocols. Arc 
different quantum protocols helped, hurt, or left intact 
by ESD? Previous results suggest a possible connection 
between the loss of certain types of entanglement in the 
four qubit cluster state and the fidelity with which mea- 
surements on the four qubit state will lead to the desired 
state on the remaining, unmeasured qubits fbl ]. The cur- 
rent paper expands these results by exploring additional 
decoherence mechanisms and calculating state indepen- 
dent accuracy measures such as the gate fidelity. Other 
explicit studies of the effect of ESD on quantum infor- 
mation protocols include the three-qubit phase flip code, 
a four qubit decoherence free subspace and a three qubit 
noiseless subsystem [TH, [H| ■ None of these studies find a 
correlation between the accuracy of the protocol imple- 
mentation and the advent of ESD. 
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II. CLUSTER STATES 



III. ENTANGLEMENT MEASURES 



The cluster state [13] is a specific type of entangled 
state that can be used as a resource for measurement- 
based quantum computation Q. A cluster state can be 
constructed by rotating all qubits into the state |+) = 
-^(|0) + |1)) and applying control phase (CZ) gates be- 
tween desired pairs. In a graphical picture of a cluster 
state, qubits are represented by circles and pairs of qubits 
that have been connected via CZ gates are connected by 
a line. A cluster state with qubits arranged in a two- 
dimensional lattice, such that each (non-edge) qubit has 
been connected via CZ gates with its four nearest neigh- 
bors, suffices for universal QC. 

After constructing the cluster state, any quantum com- 
putational algorithm can be implemented using only 
single-qubit measurements performed along axes in the 
x-y plane of the qubit, i.e. the plane spanned by |+) = 
^(|0) + | +*} = ^=(|0) + These process- 

ing measurements are performed by column, from left to 
right, until only the last column remains unmeasured. 
The last column contains the output state of the quan- 
tum algorithm which can be extracted by a final readout 
measurement. One can view each row of the cluster- 
state lattice as the evolution of a single logical qubit in 
time. Two (logical) qubit gates are performed via con- 
nections between two rows of the cluster state. CZ gates 
in particular are 'built-in' to the cluster state and simple 
measurement on two connected qubits in different rows 
automatically implements the gate. 

Measurement of a physical qubit in the cluster state at 
an angle 4> from the x-axis in the x-y plane implements 
a rotation on the logical qubit given by X(Ttm) HZ (4>), 
where H is the Hadamard gate and Z(a) (X(a)) is a z- 
(x-) rotation by an angle a [l8|. The dependence of the 
logical operation on the outcome of the measurement is 
determined by the value of to = 0,1 for measurement 
outcome — 1,+1, respectively. An arbitrary single log- 
ical qubit rotation can be implemented via three such 
measuremnts yielding: 

HZ(9i + ■Kme 1 )X(9 2 + irmg 2 )Z(9 3 + ttttiq 3 ) , 

where (6±, 82, #3) are the Euler angles of the rotation. As 
an example, by drawing the Euler angles according to 
the Haar measure, a random single-qubit rotation can be 
implemented. 

We explore an arbitrary single (logical) qubit cluster- 
based rotation performed on an arbitrary initial state 
in a decohering environment. To construct the rele- 
vant cluster, a qubit is placed in the desired initial state 
\ipi n (a, /3)) = cosa|0) + e 4/3 sina|l), where p,„(a, (3) = 
\ipin{ct, /3))(ipi n (a, /3)\. Three additional qubits (num- 
bered 2-4) are rotated into the |+)-state and CZ gates 
are then applied between the original qubit and 2, 2 and 
3, and 3 and 4. The four qubit initial (pure) state is 
thus \^4i(a,j3)) = CZ 3A CZ 23 CZ l2 (\^ in (a,fi)) ® |+)® 3 ) 
or p±i{a,l3) = |-0 4 /(a,/3))(-04/(a,/3)|- 



To quantify and monitor entanglement in the above 
constructed types of cluster states as they undergo deco- 
herence we use an entanglement measure known as the 
negativity, N , defined as the most negative eigenvalue of 
the parital transpose of the system density matrix 
There are a number of inequivalcnt forms of the nega- 
tivity for any four qubit system: the partial transpose 
may be taken with respect to any single qubit, or 
the partial transpose may be taken with respect to any 
two qubits: N^' k \ The negativity thus defined does 
not differentiate different types of entanglement. Fur- 
thermore, due to the possible presence of bound entan- 
glement, the disappearance of all measureable negativity 
does not guarantee that the state is separable. However, 
the presence of negativity does ensure the presence of 
distillable entanglement in the system. 

A method of monitoring specifically four qubit clus- 
ter type entanglement is via the expectation value of the 
state with respect to an appropriate entanglement wit- 
ness poj . Entanglement witnesses are observables with 
positive or zero expectation value for all states not in a 
specified entanglement class and a negative expectation 
value for at least one state of the specified entanglement 
class. Entanglement witnesses may allow for an efficient, 
though imperfect, means of experimentally determining 
whether entanglement is present in a state (as opposed to 
inefficient state tomography). This is especially impor- 
tant for experiments with any more than a few qubits as 
it may be the only practical means of deciding whether 
or not sufficient entanglement is present in the system. 
The entanglement witnesses used here are specifically de- 
signed to detect four qubit cluster type entanglement of 
the kind exhibited by states of the form p4 L i{a 1 (3). In 
Ref. 21[ an entanglement witness is constructed for a 
cluster state in which the first qubit is |+). It is given 
by W+ = 1/2 — (04/(71-/4,0). For the current study we 
modify this witness by a phase rotation of angle f3 on the 
first qubit yielding witnesses of the form: 



Wp = 1/2 - e.-' i ' 3 ^ /2 j O4/(7r/4,0)e 1 ' 3 ^/ 2 , 



(1) 



where crj is the Pauli z spin operator on qubit k and (3 
is the phase of the initial state \ipi n {a, /?)). This witness 
more accurately determines whether the cluster states of 
interest in this work have any four-qubit cluster entan- 
glement. 



IV. SUPEROPERATOR REPRESENTATION 

We would like to completely describe the evolution of 
the single logical qubit undergoing an arbitrary cluster- 
based rotation in the presence of decoherence. To do 
so we need to account for both the decoherence and the 
measurements of the (three) physical qubits. For this 
study we assume that there is no interaction between 
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the qubits of the cluster state (beyond the initial condi- 
tional phase gates used to construct the cluster state). 
We further assume that all decoherence occurs after con- 
struction of the cluster state but before measurements. 
Measurement is done on each of the first three qubits in 
bases at angle Oi, i = 1,2, 3, from the positive x-axis. As 
noted above, the measurement bases are chosen so as to 
implement the desired logical qubit rotation. After mea- 
surement the final state of the logical qubit resides on the 
fourth, unmeasured, physical qubit and is a function of 
the initial state of the logical qubit \i(j in (a, ft)), the deco- 
herence strength, p, and the three measurement angles, 
Oi- Pout = Pout (a, ft,p, 0i, 2 , 3 ) . 

To construct the dynamical superoperator of the one 
qubit logical gate we follow the method described in 
[H,[22j]. We construct the appropriate cluster, apply deco- 
hering evolution, and perform the desired measurements 
on a set of states, ^^(a, ft)), which span the one qubit 
Hilbert space (Hilbert space dimension TV = 2). From 
this we can construct the N 2 x N 2 Liouvillian superop- 
erator, S, where 

S(p, 0i, 2 , 3 )p m (a, ft) = p out (a, ft,p, 0i, 2 , 3 ). (2) 

Note that in Liouvillian space, density matrices are col- 
umn vectors of dimension N 2 x 1. From the superopera- 
tor S we can construct the corresponding N x N Kraus 
operators following [23[ . An analysis of the Kraus opera- 
tor representation of the scenarios outlined below is done 
in the Appendix B. 

A. Accuracy Measures 

There are two accuracy measures that we find useful 
for our analysis and that we use to compare the accuracy 
of the implemented gate to the evolution of the entan- 
glement. These measures quantify how well the system 
performed the desired operation and are thus vital in 
experimental work. The first accuracy measure we uti- 
lize is the cluster state fidelity of the four qubit state, 
P4i?(a, ft,p), before measurement but after decoherence 
as a function of p. This is given by: 

F c = Tr\p iI {a,ft)p iF {a,ft,p)% (3) 

This is a simple measure which tells how close the ac- 
tual final state is to the desired one in the presence of 
decoherence. The second accuracy measure is the gate 
fidelity of the attempted single (logical) qubit rotation, 
C/(0i,02,03). The gate fidelity quantifies the accuracy 
with which the attempted evolution was achieved inde- 
pendent of the initial state of the system. The superop- 
erator allows us to calculate the gate fidelity via: 

F<> = Tr[S(0, 0i, 2j e 3 )S(p, 0i, 2 , 3 ) f ]. (4) 

where 

5(0, 0i, 02, 3 ) - U{Bi, 2 , 3 ) ® Conj(U(0i, 2 , 3 )). (5) 



In the next three sections we look at decohering envi- 
ronments of experimental interest: phase damping, am- 
plitude damping and depolarization. In all three we ex- 
plore the entanglement evolution as a function of deco- 
herence strength assuming that the decoherence occurs 
prior to measurement. Measurements are performed on 
the decohered state and thus the evolution of the fidelity 
of the implemented operation can be compared to the 
evolution of the entanglement. Our goal is to see what 
correlations exist between entanglement degradation and 
the accuracy with which the cluster state can be used to 
implement the desired single logical qubit rotation. We 
will also note occurrences of ESD and what effect this 
phenomenon may have on the ability of the system to 
implement the desired rotation. 

V. DEPHASING 

The first decohering environment we explore is inde- 
pendent qubit phase damping. This environment can be 
completely defined by the Kraus operators 

where we have introduced the dephasing parameter p. 
When all four qubits undergo dephasing we have 16 
Kraus operators each of the form A m = {Ki <E) Kj <£> 
Kk <8> Kg) where m = 1,2,..., 16 and i,j,k,£ = 1,2. 
The final state after dephasing is thus, p^p{a, ft,p) = 
Sm=i A m p±i{oL, ft)A\ n . All of the below calculations (for 
this and the other decohering environments) are done 
with respect to p, where the exact behavior of p as a 
function of time is left unspecified so as to accomodate 
various possible experimentally relevant behaviors. As 
an example, one might have p = 1 — eT KT where r is time 
and k is the decay constant. In the case of independent 
qubit dephasing, this decoherence behavior would decay 
off-diagonal terms of the density matrix as a power of 
e~ KT and thus go to zero (i. e. p — > 1) only in the limit 
of infinite times. We also assume equal time dephasing 
for all four qubits. 

The evolution under dephasing of entanglement in the 
four-qubit cluster state, p^i, is a function of the dephas- 
ing strength, p, and initial state parameterzied by a, 
while we find that ft has practically no effect on the en- 
tanglement. Figure Q] illustrates the evolution of vari- 
ous negativity measures and the expectation value of the 
entanglement witness. The system exhibits ESD with 
respect to and M 1 * 2 ) at p = 2(^2 - 1) ~ .828, 

and ESD with respect to other negativities at p ~ .938. 
Four-qubit cluster entanglement detected via the entan- 
glement witness, W/3, disappears much more quickly, at 
decoherence strengths p < .5. 

We numerically determine, S z , the superoperator de- 
scribing the evolution of the single logical qubit in a de- 
phasing environment, via Eq. [5] and show the results in 
Appendix A. S z provides a complete description of the 
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FIG. 1: Negativity as a function of dephasing strength, p, and 
initial state (paramaterized by a; ft has little effect on any 
of the entanglement measures and thus we set ,8 = 0): (a) 
N"', partial transpose taken with respect to the first qubit, 
(b) N^ 1 ' 2 ' , partial transpose taken with respect to the first 
two qubits, (c) N*- 1 ' 3 ' partial transpose taken with respect to 
qubits one and three, (d) Evolution of the expectation value 
of the entanglement witness as a function of initial state (the 
expectation value is not dependent on /?) and decoherence 
strength. Notice that the dephasing strength at which the 
expectation value goes to zero is dependent on a and is well 
below the point where ESD is exhibited for iV*- 1 ^. The four 
qubit cluster entanglement can only be observed at low levels 
of decoherence, p < .5. 

single logical qubit evolution including the effects of the 
independent qubit dephasing environment and the mea- 
surements necessary to implement the logical qubit rota- 
tion. From S z we can determine the output state for any 
input state and we can compute the gate fidelity of an 
arbitrary cluster-based single qubit rotation as a function 
of dephasing strength: 

F% = — (lO + 6p + p(p-6p-7) + q(p + 2p-2)cos20 2 
16 

+ 2q(p + 2p- 2)cos6»2Cos26» 3 ), (7) 

where p = ^/l — p. We can also determine the fidelity 
measure between the decohered four qubit cluster state, 
and the initial cluster state. This fidelity measure 
quantifies how much the initial state, before measure- 
ment, is degraded by decoherence, and, as a function of 
dephasing strength is given by: 

Fz = ^r(16(l+p)+p(j5-6p-14)-p(p-2p-2)cos4a). 

(8) 

These fidelity measures arc plotted in Fig. [3J 
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6, a 

FIG. 2: (a) Gate fidelity of an arbitrary single qubit rota- 
tion as a function of the rotation (paramaterized by Euler 
angle 82 with #3 set to zero; note that the gate fidelity is in- 
dependent of 81), and dephasing strength p. (b) Fidelity of 
pre-measurement four qubit cluster state as a function of ini- 
tial state paramaterized by a (this fidelity is independent of 
ft) and p. 

From the figures we see a slight similarity in behav- 
ior between A' 1,2 ) and the four-qubit state fidelity. Be- 
yond that there is no correlation between the evolution of 
any of the other entanglement measures and the accuracy 
with which the logical qubit rotation is performed or the 
fidelity of the decohered four qubit state. Importantly, 
we see no manifestation of the ESD in the accuracy mea- 
sures. There is no sudden change of behavior or sharp 
decline at the decoherence strength where ESD occurs. 



VI. AMPLITUDE DAMPING 

We now turn to an independent qubit amplitude damp- 
ing environment and, as above, explore the entanglement 
evolution and compare it to various accuracy measures. 
The Kraus operators for this environment are: 

*-(i^).*-(2tf). (9) 

where the amplitude damping strength is again denoted 
V- 

The evolution of various negativity measures is shown 
in Fig. [3] As in the case of the dephasing environment 
the variable (3 has no apparent effect on the negativity. 
We first note that, unlike the dephasing environment, en- 
tanglement sudden death (ESD) does not occur for any 
entanglement measure. The negativity may reach zero 
only in the limit of p approaching 1. However, the en- 
tanglement witness detecting four qubit cluster entangle- 
ment goes to zero rather quickly, a maximum amplitude 
damping strength of p = .2 for (3 = and increasing to 
at most p ~ .3 at /? ~ 7r/3. This shows a quick demise 
specifically for the four-qubit cluster entanglement (or 
indicates the inefficiency of the witnesses). 

We determine Sa, the superoperator of an arbitrary 
single qubit rotation in the presence of amplitude damp- 
ing, via Eq. [2] and show the results in Appendix A. Sa 
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FIG. 3: Evolution of various negativity measures as a func- 
tion of initial state (parameterized by a; /3 has little effect 
on any of the entanglement measures and is set to zero) and 
amplitude damping strength, p. (a) N^ 1 ', the entanglement 
goes to zero at a = tt/2. (b) iV (1 ' 2) , (c) iV (1 ' 3) , for these 
measures the entanglement goes to zero only in the limit of 
p — > 1. (d) Expectation value of four qubit cluster state with 
respect to entanglement witness W/3, with /3 = 0, as a func- 
tion of initial state and decoherence strength. The four qubit 
cluster entanglement goes undetected at very low decoher- 
ence strengths (p < .2) despite the presence of some sort of 
entanglement for any non-zero p. 

gives a complete description of the single logical qubit 
evolution including the effects of the independent qubit 
amplitude damping environment and the measurements 
necessary to implement the logical qubit rotation. Sa al- 
lows the determination of the output state for any input 
state and can be used to calculate the gate fidelity with 
which any desired rotation can be impelemented. This is 
given by: 

F 9 A = i(10 + 6p + p(p-6p-7) + g(p + 2p-2)cos26» 2 

+ 2q(p + 2p- 2)cos6>|cos26> 3 ). (10) 

This is in fact the identical gate fidelity calculated for the 
dephasing environment, = F-?. This is not surpris- 
ing considering that the Kraus operators governing the 
decoherent environments are very similar. Despite hav- 
ing the exact same gate fidelity, the cluster state fidelity 
measure between the decohered four qubit cluster state, 
P4F, and the initial state (which is a function of a and 
p only) in the amplitude damping environment is not at 
all similar to that of the dephasing environment, exhibit- 
ing a region around a = ir/2 where the fidelity attains 
extremely low values. This demonstrates the importance 
of utilizing multiple accuracy measures. Both of these 
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FIG. 4: (a) Gate fidelity as a function of amplitude damp- 
ing strength and choice of rotation (paramaterized by 82 with 
83 — 0; note that the gate fidelity is independent of Si), (b) 
Contour plot of four-qubit state fidelity as a function of am- 
plitude damping strength and initial state (paramaterized by 
a; the fidelity is independent of /3). 

meaures are shown in Fig. |4j 

Though the gate fidelity of an arbitrary single qubit 
cluster-based rotation is the same for both a phase 
damping environment and an amplitude damping envi- 
ronment, the amplitude damping environment does not 
cause ESD in any entanglement metrics. This lack of 
ESD suggests that dephasing is more harmful to the en- 
tanglement types found in the four qubit cluster state 
than is amplitude damping. Thus, rotating the phys- 
ical system qubits, such that a dephasing environment 
acts as an amplitude damping environment could con- 
serve entanglement despite not increasing the accuracy 
of the implementation of the single logical qubit rota- 
tion. Furthermore, the lack of ESD in only the amplitude 
damping environment highlights the difference in entan- 
glement evolution in the two environments (though the 
behavior is similar) again demonstrating the lack 
of correlation between entanglement evolution and the 
accuracy of the implemented protocol. 

VII. DEPOLARIZATION 

The final decohering environment in our study is 
the independent qubit depolarizing environment and, as 
above, we explore the entanglement evolution and com- 
pare it to various accuracy measures. The Kraus opera- 
tors for this environment are: 

V ° V^f/ 2 

where the <jg are the Pauli spin operators, £ = x,y, z and 
p is now the depolarizing strength. 

The evolution of entanglement in the four-qubit clus- 
ter state, pn, in a depolarizing environment depends on 
the depolarizing strength, p, and initial the state a (the 
entanglement appears to be independent of j3). Figure 
[5] illustrates the evolution of various negativity measures 
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FIG. 5: Negativity as a function of depolarizing strength, p, 
and initial state paramaterized by a: (a) N^ 1 ', partial trans- 
pose taken with respect to the first qubit, (b) iV*- 1 ' 2 ', partial 
transpose taken with respect to the first two qubits, (c) jV"' 1 ' 3 ) 
partial transpose taken with respect to qubits one and three, 
(d) Evolution of expectation value of entanglement witness as 
a function of initial state and decoherence strength. Notice 
that the evolution is similar to iV' 1 ) though the expectation 
value goes to zero well before ESD of iV' 1 ' The four qubit 
cluster entanglement can only be observed at low levels of 
decoherence, p < .2. 
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FIG. 6: (a) Gate fidelity for arbitrary single qubit rotation 
in depolarization environment as a function of depolarization 
strength and the rotation (paramaterized by 62, in this case 
the gate fidelity does not depend either on 9\ or 63). (b) 
Fidelity of four qubit cluster state as a function of initial 
state (parameterized by a as Fp does not depend on /3) and 
depolarizing strength. 

Both of these meaures are shown in Fig. [6] and decrease 
much more quickly in a depolarizing environment than 
in dephasing or amplitude damping environments. 

Once again we find only superficial correlations be- 
tween the entanglement evolution and either the gate 
fidelity or the fidelity of the cluster state prior to mea- 
surement. These superficial correlations do not give rise 
to any problems regarding the viability of quantum com- 
puting. While clearly both entanglement and fidelity de- 
crease as the decoherence strength increases, we do not 
find any signature of ESD in the fidelity functions. 



and the expectation value of the entanglement witness. 
All of the negativity measures exhibit ESD at p < .45. 
The expectation value of Wp which enables the detection 
of four-qubit cluster entanglement, disappears at deco- 
herence strengths p < .2, though its general behavior 
appears similar to that of . 

We determine Sp, the superoperator of the arbitrary 
single logical qubit cluster based rotation in a depolariz- 
ing environment, via Eq. [5] and show the results in Ap- 
pendix A. This gives a complete description of the single 
logical qubit evolution including the effects of the inde- 
pendent qubit depolarizing environment and the mea- 
surements necessary to implement the logical qubit ro- 
tation. Sp allows the determination of the output state 
for any input state and can be used to calculate the gate 
fidelity of an arbitrary single qubit rotation: 

Ff, = I((p-2)(-4+p(7+p(p-6)))+9y cos20 a ). (12) 

The fidelity measure between the decohered four qubit 
cluster state, p4p, and the initial state is a function of a 
and p only and is given by 

Fp = ^(P - 2) 2 (3p(3p - 5) + 8 - qpcosAa). (13) 



VIII. CONCLUSION 

In conclusion, we have studied the evolution of entan- 
glement in a four-qubit cluster state as a function of ini- 
tial state and decoherence strength for three different de- 
cohering environments: dephasing, amplitude damping, 
and depolarization. This behavior was compared to the 
fidelity with which the cluster-based arbitrary one-qubit 
rotation could be performed on the decohered qubits. In 
addition, we commented on occurrences of entanglement 
sudden death and noted that this phenomenon is not rel- 
evant to the fidelity with which the intended gate is im- 
plemented. This point should be noted when discucssing 
the role of entanglement in quantum computation. We 
believe that studies such as the one presented here are 
essential for planning experimental implementations of 
cluster based quantum computation and for calculation 
of cluster based fault tolerance thresholds. 
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Appendix A: Arbitrary Single Qubit Rotation 
Superoperators 

In this Appendix we provide the expressions for super- 
operators describing the evolution of a single logical qubit 



in a cluster based quantum computation attempting to 
implement an arbitrary rotation described by Euler an- 
gles (9i,8 2 ,0 3 ) in the different decohering environments. 
The superoperator for the case of qubits in a dephasing 
environment is given by 



/ (1 - qs2s3) 



l q(c3 — ipc2s3) 



q(c2 - ipc3s2) e ll>1 q(qc2c3 + ip(s2 + s3)) -i 
q(c2 + ipc3s2) -e l8l q{qc2c3 - ip(s2 - s3)) e 
V (1 + qs2s3) e l6l q(c3 - ipc2s3) 



-e- Wl q(c3 + ipc2s3) 
Wl q{qc2c3 + ip(s2~- s3)) 
h q(qc2c3 - ip(s2 + s3)) 
e- iBl q(c3 + ipc2s3) 



(1 + qs2s3) 
—q{c2 — ipc3s2) 
—q(c2 + ipc3s2) 

(1 - qs2s3) 



(Al) 



where q = p — 1 , p = y/1 — p, and we write sin 9 j and cos 9 j for j = 1 , 2 , 3 as sj and cj . 
The superoperator for the amplitude damping environment is given by 



S 



( (l+p) +q 2 s2s3 
q(c2 — ipc3s2) 
q{c2 + ipc3s2) 

\ -q(l+qs2s3) 



e t8l q 2 (c3-ipc2s3) 
e Wl q(qc2c3 + ip(s2 + s3)) 
-e l6l q{qc2ci - ip(s2 - s3)) 
-e iei g 2 (c3-ipc2s3) 



2 



q 2 {c3 + ipc2s3) (1 + p) - q 2 s2s3 \ 
9l q{qc2c3 + ip(s2 - s3)) -q(c2 - ipc3s2) 
l q(qc2c3 - ip(s2 + s3)) -g(c2 + ipc3s2) 



il q 2 (c3 + ipc2s3) 



(A2) 



-g(l - gs2s3) J 



We note that the second and third rows of the amplitude damping superoperator are exactly the same as the second 
and third rows of the dephasing superoperator. 

The superoperator for the depolarizing environment is given by 



/ 1 - q 3 s2s3 



-e- l01 q 3 {c3 + iqc2s3) 



S P = ^ 



- l8l q 3 (-c3 + iqc2s3) l + q 3 s2s3 

Wl q 3 {qc2c3 + i{s2 - s3)) q 2 {c2 + iqc3s2) 
-q 2 (c2-iqc3s2) -e Wl q 3 (qc2c3 - i(s2 - s3)) e- i6l q 3 (qc2c3 - i(s2 + s3)) g 2 (c2 - iqc3s2) 



-q 2 (c2 + iqc3s2) e i9l q 3 {qc2c3 + i(s2 + s3)) -e 



\ 1 + g J s2s3 



e- iei g J (c3 + ^c2s3) 



-ifli„3 



g J (c3 — iqc2s3) 



(A3) 



1 - <fs2s3 / 



The above superoperators promise to be useful for experimental realizations of this cluster state protocol as they 
can be used to characterize a given environment. 



Appendix B: Kraus Operator Representation 



complete evolution of the system is given by 



In the main part of this paper we determined the su- 
peroperators for single logical qubit rotations in a cluster 
based quantum computer undergoing different types of 
decoherence. Recasting these superoperators in terms of 
Kraus operators gives additional insight into the evolu- 
tion of the logical information under the arbitrary qubit 
rotation as a function of decoherence. To calculate the 
Kraus operators from the superoperator one first deter- 
mines the Choi matrix. Each Kraus operator K a is a 
Choi matrix eigenvector (unstacked so that its dimen- 
sion is iV x N), times the square root of the correspond- 
ing Choi matrix eigenvalue divided by ./V [23[ ■ We define 
the amplitude of a given Kraus operator, A ai to be the 
square root of the Choi matrix eigenvalue divided by N, 
A a = \J \ a /N. The higher the amplitude of a Kraus 
operator the more significant its effect on the overall sys- 
tem dynamics. This method of Kraus operator construc- 
tion maximizes the amplitude of one (and hence the most 
significant) Kraus operator. Using Kraus operators, the 



Y, K a (j>, 9 ll 9 2l 9 3 ) Pm {a, P)Ka{p, 0^62,6^ = Pout . 

a 

(Bl) 

Clearly, if there is only one Kraus operator it will be 
unitary with an amplitude of 1. In this way, unitarity 
of the evolution can be quantified by A\ , the amplitude 
of the first Kraus oprerator. In addition, the accuracy 
of the applied evolution can be quantified by the fidelity 
or correlation of the first Kraus operator as compared to 
the desired unitary [22|]. The fidelity is given by: 



and the correlation is given by: 



C 1 = 



Tr[j7t^i 



Tr[C/t[/]Tr[A' 1 t A'i] 



(B2) 



(B3) 



The fidelity measure accounts for both decoherent losses, 
a change in the purity of the state, and coherent errors, 
what we might call a change in the 'direction' of the state, 



8 




FIG. 7: (left) Kraus operator amplitudes as a function of de- 
phasing strength. The highest Kraus operator amplitude de- 
creases linearly until becoming equal to the amplitudes of the 
other three Kraus operators. The behavior of the highest and 
lowest amplitudes are independent of all measurement angles 
while the middle two amplitudes depend slightly on 62 and 
6*3. (right) Fidelity (solid line) and correlation (dashed line) 
of first Kraus operator as a function of dephasing strength. 




FIG. 8: (left) Kraus operator amplitudes as a function of 
amplitude damping strength. The behavior of all amplitudes 
depend slightly on 62 and 83. (right) Fidelity (solid line) and 
correlation (dashed line) of first Kraus operator as a function 
of dephasing strength. Note that in this case the correlation 
does not remain constant implying a coherent affect on the 
dynamics of the system due to amplitude damping. 

effected by the first Kraus operator. The correlation is 
unaffected by a change in magnitude of the state. 

We start with the Kraus operators of the dephasing 
environment. The amplitude of the Kraus operators as a 
function of decoherence strength is shown in Fig. [7] The 
amplitude of the first Kraus operator decreases linearly 
and the amplitude of the other three Kraus operators 
increase until, at p = I, the four Kraus operators have 
equal amplitudes. At that limit each of the four Kraus 
operator matrices have an element equal to 1/V% in one 
of the corners and all the other elements are zero. The 
fidelity of the first Kraus operator also decreases linearly 
with dephasing strength while the correlation remains 
constant at I (until very high p) implying purely deco- 
herent evolution. 

We noted in the main part of the paper that the gate 
fidelities of a single logical qubit cluster state-based arbi- 
trary rotation in a dephasing environment and amplitude 
damping environment are the same. Nevertheless, we 



find that their Kraus operator representations are very 
different. Two of the Kraus operator amplitudes in an 
amplitude damping environment go to zero as p — > 1. 
The remaining two Kraus operator matrices have a one 
in the upper right or left corner and zeros elsewhere. For 
values of p < 1 , the amplitude of the first Kraus oper- 
ator decreases faster in the amplitude damping environ- 




FIG. 9: (left) Kraus operator amplitudes as a function of 
depolarizing strength. The behavior of all amplitudes depend 
slightly on 62- (right) Fidelity (solid line) and correlation 
(dashed line) of first Kraus operator as a function of dephasing 
strength. 



ment than in the dephasing environment. However, this 
descent slows as p approaches one. The second Kraus 
operator always plays a more significant role in the am- 
plitude damping environment than in the dephasing en- 
vironment. 

The fidelity of the first Kraus operator as function of 
amplitude damping strength p decreases linearly (faster 
than the dephasing environment) before rounding off at 
high values of p while the correlation decreases, staying 
near one only at low values of p. This behavior, portrayed 
in Fig.0 suggests that amplitude damping, despite being 
decohcrent dynamics, has a coherent affect on the system 
dynamics. Rotating the system so that the amplitude 
damping acts as phase damping may increase the fidelity 
and correlation of the first Kraus operator but will not 
increase the gate fidelity of the attempted logical qubit 
rotation. 

For all the decohering environments studied here, the 
first Kraus operator amplitude decreases fastest (and not 
linearly) in a depolarizing environment. The increase of 
the lowest amplitude Kraus operator is also not linear. 
However, in the limit of p — > 1, the depolarizing envi- 
ronment is like the dephasing environment in that the 
amplitudes of all four Kraus operators converge to .5, as 
demonstrated in Fig. |H1 The fidelity in a dcpolarzing en- 
vironment also decreases faster than the other decoherent 
environments while the correlation remains constant at 
one, demonstrating that the evolution is entirely deco- 
herent. 
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